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Abstract 

A neural network method is developed to discriminate direct photons from the neutral 
pion background in the PHOS spectrometer of the ALICE experiment at the LHC 
collider. The neural net has been trained to distinguish different classes of events by 
analyzing the energy-profile tensor of a cluster in its eigen vector coordinate system. 
Monte-Carlo simulations show that this method diminishes by an order of magnitude 
the probability of 7r°-meson misidentification as a photon with respect to the direct 
photon identification efficiency in the energy range up to 120 GeV. 



1 Introduction 

The Large Ion Collider Experiment ALICE [1] is intended to study heavy-ion col- 
lisions at the energy of 5.5 TeV per nucleon at the LHC collider. For photons detection 
ALICE is equipped by the PHOton Spectrometer (PHOS) [2] which is a high-granularity 
electromagnetic calorimeter built of lead-tungstate crystals (PbWOzt) with the transverse 
dimension of 2.2 x 2.2 cm 2 and the length of 18 cm. The PHOS detector consists of five 
modules, each made as a matrix of 56 x 64 crystals located at 460 cm from the beams 
interaction point. The spectrometer is positioned at the bottom of the ALICE set-up 
covering \r)\ < 0.135 in pseudorapidity and 100° in the azimuthal angle. 

The ALICE set-up is rather transparent to 7-radiation. To reach PHOS, the pro- 
duced photons pass through the Inner Tracking System (ITS) [3], the Time Projection 
Chamber (TPC) [4] and the Charge Particle Veto Detector (CPV) [5]. The stratum of 
medium takes ~ 0.1 of the radiation length [6] (with the main contribution from ITS and 
TPC), which gives a small probability of the secondary photon production in the medium. 

One of the tasks of the ALICE experiment is to detect direct photons carrying the 
information about fundamental processes at extreme conditions of the quark matter. The 
energy range of such photons extends from 0.4 GeV to 10 GeV for thermal radiation of 
Quark Gluon Plasma (QGP) and higher, up to hundreds GeV, for radiation occurring also 
in early collisions of the constituents. In the latter case the essential background is arising 
from the two-photon decay of 7r°-mesons, produced at the same energies as photons, due 
to the merge of decay photons into one shower in PHOS. Contribution of the photonic 
decay from heavier meson resonances (p, rj, u, etc. ) exists as well. Such photons provide 
a rather heavy background to direct photon production over the whole Pt range and they 
should be subtracted from the data. It leads to the necessity of identifying the PHOS 
showers as those produced by the photons or 7r°-mesons. 

The straightforward way to discriminate the considered background from direct 
photons is to exploit a powerful reconstruction program tuned especially for the given 
PHOS structure and able to resolve the fine problem of the "direct" photon selection 
(say, by the maximum likelihood method), and correspondingly, to find the real number 
of the showers even in the overlapped clusters. To accomplish this task, a program should 
correctly take into account the shower profiles of photons and charged particles at dif- 
ferent incident angles, as well as their fluctuations, the electronics noise, the threshold 
on registered signals (that increases the number of zero channels in the selected cluster), 
etc. In practice such a perfect program is often too complicated for fast realization with 
long tuning and commissioning and additionally needs essential computational resources. 
Therefore, real reconstruction programs are usually created in some simplification as- 
sumptions which results in the degradation of its discrimination power for direct photon 
selection. Such simplified programs are often based on the recognition of the shapes of 
the showers produced by different kinds of particles. 

In this paper we apply the Artificial Neural-Network (ANN) approach [7] for the 
direct photon identification in PHOS. The ANN-method has already recommended itself 
as a powerful tool in different applications of high energy physics, e.g. quark and gluon jet 
separation, b-quark identification, Higgs boson search, selection of the rare neutral meson 
decay modes violating C-parity, etc. [8-13]. 

The main peculiarity of our method is the use of the energy-profile tensor of a 
cluster, which components are calculated in its eigen-vector coordinate system, with the 
aims of the neural net training and the afterward event classification. Our calculations 
show essential recognition capacities of this procedure that were examined using a sample 



of Monte-Carlo generated events simulating the isolating production of direct photons 
and 7r°-mesons for the real ALICE set-up. 

2 Application of neural network algorithm method 

In the analysis of experimental data, a standard procedure of selecting signal events 
is based on various cuts of observed kinematics variables. A general case of such cuts 
corresponds to a particular set of functions, called feature functions or neurons. In general 
terms, neural networks are represented as a large number of interconnected functional 
units named nodes with a neuron in each of them. The data processing is organized in 
the most common, to-date, architecture called the Multilayer Perception (MLP). MLP 
incorporates one input layer through which the initial data (features) are injected, several 
hidden layers, and one output layer of neurons. 

Output response Oi of the i-th neuron is obtained by the so called sigmoid function 
/ (0 < / < 1) dependent on the weighted sum of all input signals Sj to this neuron: 

o i = f(£w ir s j + e i ) , (l) 
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where index j runs over all the inputs of the i-th node, Wij and 6i are the weights and the 
correspondent thresholds optimizing the selecting power of the procedure. Function / is 
defined as 

f( x ) = (l+tanh(x))/2 . (2) 

The output layer, consisting of one node, provides a neural-net response Snn, wedged 
between and 1, and used to classify the events. 

The use of the neural network is a two-step process, i.e. a learning stage followed 
by an application stage. During the learning phase using the Monte-Carlo simulation, we 
know about every event whether it is a background or a signal one. The optimal values 
of weights Wij and thresholds (see (1)), for resolving the problem are determined by 
minimizing the functional L(wij,6i) 

1 N (k) (k) 
L(Wij, ^i) = — XI l^lcarn ~ $Nn\ j (3) 
Z k=l 

where index k runs over all N training events, S^ rn = for the background and S^ rn = 1 
for desired events, S$ N are the ANN response. The details of the whole minimization 
procedure can be found in [14]. And finally the quality of learning is tested with an 
independent sample of Monte-Carlo events. 

2.1 Input variables for the ANN 

Reconstruction programs for cellular electromagnetic calorimeters provide data as a 
set of clusters defined as a group of cells with a common edge. Every cluster is characterized 
by the amplitudes of signals from the cells and coordinates of the cluster cells on the 
detector plane. The total length of the data array for one cluster is 3N C , where N c is a 
number of cells in the cluster. This array contains exhausting experimental information 
about the cluster but, however, it is not convenient for the aim of the ANN approach, at 
least due to the varying data length from event to event. To overcome the latter problem 
one can use only the limited number of major cluster cells (say, 3x3 around the cell with 
the maximal amplitude [15] ), but it definitely leads to the uncontrolled information loss. 



The essential part of our ANN- approach to 7/71° selection is to choose such a fixed 
number of variables that carry, nevertheless, the principal part of information with its 
volume big enough to find a solution and which length is independent of the cluster size. 
First of all each cluster is characterized by a position of its center-of-gravity X g as well as 
by the center of the cell with the maximal signal in the cluster X m . Vector X g is defined 
according to the standard formula 

X g = E*" )E k/E E k> ( 4 ) 
k k 

where index k runs over all cluster cells, x[ k ^ and being the coordinates of the fc-th 
cluster cell in an arbitrary coordinate system. 

Then we introduce an energy-profile tensor Qij of a cluster (that can be also in- 
terpreted as a 2 x 2 matrix) calculated in the local cluster coordinate system. There are 
two natural possibilities to set the origin of this coordinate system: either in the center of 

— * — * 

gravity X g or in the center of the cell with the maximal signal in the cluster X m . In our 
calculations we choose the latter option. After that is defined as 

Qij = Y,4 k) x? ) E k ; i,j = 1,2. (5) 
k 

Here index k runs over all the cluster cells, and x^ are coordinates of the fc-th cluster 
cell, for definition, in non-dimensional units of the crystal transverse size, and E k is the 
energy detected in this cell. Tensor (5) is a quadratic positive-definite form. 

The introduced tensor (5) reflects the cluster energy profile which depends on the 
number of showers overlapped in one cluster, as well as on inclination of the photon inci- 
dence on the detector. To avoid the latter effect or, at any rate, to decrease it we have made 
(before calculating matrix Q^) a compression of the cluster space relative to the origin 
of the local cluster coordinate system X = (x w ,X2o) along the vector directed from the 
geometrical center of the PHOS module to the point Xq. The compression coefficient has 
been defined as c = cos0, where © is the angle between the perpendicular to the PHOS 
module surface and the photon propagation direction. The coordinate transformation due 
to the compression operation has been made according to the formula 

x -> i?(-0 o ) ■ E C (Q) ■ R((f> ) ■ (x - X ) + X , (6) 

where x = (xi,x 2 ) is the vector of coordinates of an arbitrary transfered point, R and E c 
are the matrices of rotation and compression, respectively: 

/ cos0o sin O \ /cosG 0\ 

y ^ ' \ - sin 0o cos O / \ 1 J v ' 

and 0o is the polar angle of the point X in the polar coordinate system with its origin 
in the geometrical center of the PHOS module. 

One can calculate two eigen values Ai and A2 (ordered as Ai > A2) of Qij and find the 
correspondent eigen vectors e\ and e*2 (normalized to a unit), defining a new coordinate 
system, where Q^ is reduced to the diagonal form. In this new system we also define 
moments M mn : 

M mn = 5>'i fc) r • (4 fe) )" E k ; m,n = 0,1,2,..., (8) 
k 



where index k runs over the cluster cells, x l and x 2 are the coordinates of the k-th 
cluster cell in the new coordinate system. Note that M 20 = Ai, M 02 = A 2 , and M 00 = E = 
J2k Ek is the total cluster energy. Such important magnitudes as distance d between hits 
of two glued photons and their effective mass M 77 can also be expressed through M mn 
(see addendum). 

We would also like to point out an essential remark. Initially the experimental in- 
formation was carried by signal amplitudes and cell coordinates, and now it is represented 
by the introduced moments M mn . We shall construct from M mn the input vector p( m ) of 
event features for ANN. One additional angle variable, not directly expressed through the 
M mn , can be added. This is angle <fi between the eigen vector e\ and vector X directed 
from the geometrical center of the PHOS module to the cluster center 



The use of this angle, together with the coordinates of the cluster center on the detector 
plane, significantly improves the quality of 7/7r°-selection mainly in the peripheral region 
of the PHOS modules. 

We have found that there are different sets of variables that allow to construct 
effective event feature vectors. They include the total detected cluster energy, the eigen 
values Ai and A2, the moments M mn , the estimation of the effective mass M 77 and distance 
d between 77-hits expressed through M mn (see formulas (A. 4) and (A. 5) ), and the angle 
variable <fi. Coordinates of the cluster center X relative to the center of the PHOS module 
can also be added, which increases the selection power of the method. 

2.2 Strategy of the 7/7r°-selection 

In this section the algorithm of the 7/7r°-selection is discussed in detail. First of 
all we note, that at low energies there is a background from 7r° decays producing two 
separated clusters in PHOS that can be taken into account by one of the statistical 
methods, based on a good resolution of PHOS for M 77 -effective masses in this case, which 
is demonstrated below in Fig.5b. The mentioned procedures can, for example, reject the 
7T°-contribution by calculating the masses of 77-combinations and by comparing them with 
the 7r°-mass. It is possible to reduce the number of such combinations at high occupancy of 
the detector by taking into account the characteristic angle of the decay cone for photons. 
These methods are beyond the scope of the current paper. Thus, further calculations 
were fulfilled following the assumption that this type of backgrounds had been eliminated 
correctly. 

Besides, there are two other essential sources of the background due to 7r°-decays, 
when only one cluster appears in PHOS. The first one results from the detector geometry, 
i.e. when one of the decay photons from n° escapes detection in PHOS due to the limited 
acceptance, whereas the second photon hits PHOS and generates a single shower and 
therefore a single cluster. It provides actually indistinguishable from the direct photons 
sample of the background photons in the whole photon energy range. This background 
is essential at relatively small energies of neutral pions. To suppress the relative part u 
(close to one) of such a background, we demand that the cone with a specially defined 
total angle 6{u) around the photon propagation should cross the detector plane inside of 
it. The angle 9(u) can be easily estimated from the isotropy of 7r° — > 77 decay in the 
meson rest frame after the Lorentz boost to the laboratory system 



= arccos( e\, Xq/\Xq\ ). 



(9) 
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where E is the energy of the analyzed cluster, (3 ~ 1 is the 7r°-meson velocity in units of 
the light speed. Note that minimal decay cone angle # m i n = 2arccos(/9) can be obtained 
from formula (10) at u — 0. In the case when both photons fly towards the detector 
direction we arrive at u < f3. To apply the considered cut we take u = 0.95. This cut 
actually works at small energies, and mainly at the edges of the calorimeter, while at high 
energies practically all events pass it successfully. 

The second background comes from the limited spatial resolution of PHOS, i.e. when 
the overlapping showers from two-photon decays of the high-energy pions form one cluster. 
This mechanism provides the major background for direct photons at high energies. Its 
suppression is mainly fulfilled by the ANN-algorithm. The step sequence runs as follows. 

First we apply the procedure of the local peak number determination in a cluster. 
It classifies the cluster as that with two overlapping showers when two local peaks are 
observed above the electronic noise fluctuations. Such classification of clusters is a common 
algorithm in reconstruction programs for cellular detectors. 

To decrease the background from the decays 7r° — > 77, when only one photon hits 
PHOS due to the limited aperture, we apply cut (10) with the use of 6{u). The remaining 
clusters with one peak passing this 2-stage preselection are analyzed by ANN. It was 
composed of three layers (see Fig. 1): input, hidden and output. The input layer consists 




nn Figure 1: Architecture of the used ANN. 



of N nodes, where N is the dimension of vector representing the event features; the 
hidden layer is built of 2iV + 1 nodes; and finally, the one-node output layer provides the 
neural-net response Snn £ (0, 1). 

The net was trained with two samples of events containing the desired signals (i.e. 
clusters from direct photons) and the background (i.e. clusters produced by the overlapped 
photons from decay 7r° — > 27). Each sample consisted of 30,000 events (clusters) left 
after preselection. Upon training the net, we tested its efficiency on another statistically 
independent signal and background samples, each consisting of 30,000 events also left 
after the preselection. 

3 Simulation of isolating photons and 7r°-mesons 

The data taken for the analysis were simulated within the ALICE framework for 
the simulation, reconstruction and data analysis, aliroot v. 3. 06. 02 [16]. Two samples 



of events were generated, one sample containing one photon and another sample having 
one 7T° per event. Photons and 7r°'s were emitted from the interaction point with the 
uniformly distributed transverse momentum in the range of < px < 120 GeV/c within 
the solid angle defined by the uniformly distributed azimuth angle 210° < <p < 330° and 
the polar angle 80° < 9 < 100°. The solid angle of the emitted particles was chosen to be 
a little larger that that of PHOS detector, to avoid various border effects. The decay of 
7r°'s was performed by the aliroot. 

The real ALICE environment was taken into account during the particle tracking 
from the interaction point to PHOS. The following detectors and infrastructure modules 
which cover the PHOS aperture were installed: PIPE, ITS, TPC, TRD, FRAME, as 
shown in Fig.2. This environment results in particle interactions with the media and the 




Figure 2: ALICE detector used for the simulation in aliroot. 



production of secondary particles which deteriorate the response of the PHOS detector. 

The response of the PHOS calorimeter was simulated by GEANT 3.21 [17], which 
was included into the aliroot package as a particle tracking and physics interaction 
tool. The showers developed by the particles passing through the calorimeter cells, gave 
the HITS which were the deposited energy of the shower in each cell. These HITS were 
digitized, i.e. the energy amplitude of the cell was measured by the two-channel ADC, 
and the integer signal of the ADC produced DIGITS. The noise with 0" no i se = 10 MeV was 
applied to the DIGITS, after which the DIGITS were passed through the 50-MeV threshold 
cut. The remaining DIGITS with a common side were grouped into clusters. Only clusters 
with the total energy greater than 500 MeV were accepted. 



The data that passed to the Neural Network analysis contained the following infor- 
mation. The event samples (photons or 7r°'s) were identified by the file name. Events had a 
header which was characterized by the event number, the energy of the produced particle 
and the number of the found clusters. The event header was followed by the cluster data 
consisting of the cluster header with the cluster number and the number of cells in the 
cluster, and the cell information containing the list of cell positions in PHOS as well as 
the cell amplitudes. 

4 Results of the 7/7r°-selection 

Analysis of the generated data shows that only 0.04% of direct photons are recon- 
structed as two-cluster events. The energies of the additional clusters are limited, they 
are practically not greater than 1.5 GeV. 

Two-dimensional plots (Fig.3) of eigen values Ai and A 2 versus energy E of 7 or 7r°, 
demonstrate the difference between the direct photons and the background that is used 
in the process of the ANN training. Fig. 4 also shows one-dimensional distributions which 
demonstrate a rather well-expressed 7/7r°-distinction. 

Fig. 5 shows comparative results of the calculations of 77-effective masses M 77 for 
overlapped and separated clusters from 7r°-decays and the same data from the sample of 
direct photons. For the overlapped clusters we use formula (A. 5) expressing M 77 through 
the moments M mn . In both cases for the overlapped and separated clusters from 7r°-decays 
there are well observed peaks at the value of 7r -mass while application of the same formula 
(A. 5) for direct photons gives a smoothly falling dependence with the growth of mass in 
the region of the 7r -meson. Smooth solid curves in the figure show the results of Gaussian 
fits of 7r°-peaks. The average value and variation for the separated clusters are equal to 
133 MeV and 6.8 MeV, respectively. 

The main results of this article concerning the quality of the ANN training for the 
7/7r°-discrimination are presented in Fig. 6 for one of the possible feature event vectors, 
selected as (E, Ai, A 2 , M 30 , M 04 , 0), where E is the measured cluster energy. The cut on 
the S NN signal equal to 0.64. This figure shows the efficiency £(7,7) of true photon 
identification as a photon, misidentification 5(7, 7r°) of 7r°-meson as a photon in the range 
of generated energies of photons and 7r°-mesons from 3 GeV to 120 GeV, and the coefficient 
of background suppression relatively to direct photons 5(7, 7r°)/e(7, 7). One can see that 
the probability £(7, 7r°) of misidentification of a neutral pion as a photon is on the level of 
a few percent in the energy range of 3 — 120 GeV with relatively high efficiency of the true 
photon identification. The rise of the n° misidentification probability at energies below 
25 GeV is caused by the tt° — * 77 decays with one photon outside PHOS, whereas the 
decrease of the true photon identification efficiency at the energies below 15 GeV is the 
result of the true photon cut decreasing the background from these decays of 7r°'s with 
only one photon inside PHOS. 

We compare our results with the data of work [18] where the coefficient 5(7, 7r°) 
was estimated for the STAR experiment as 0.15 at 20 GeV and 0.45 at 40 GeV at fixed 
£ (7>7) — 0-8- The efficiency of the 7/7r°-recognition was also calculated within the neural 
network approach for the CMS experiment in note [15]. The obtained values 5(7, ir°) varied 
from 0.25 to 0.55 at E = 20 GeV and from 0.40 to 0.55 at E = 100 GeV, dependent on 
the rapidity range, while the efficiency of the single photon recognition was kept at 91%. 
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Figure 3: Two-dimensional plots of Ai (a), (b) and A2 (c), (d) vs energy, (a) and (c) are 
direct photons, (b) and (d) are the background. 
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Figure 4: Ai- and A2 distributions for direct photons (solid lines) and background (dotted 
lines) . 
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Figure 5: a) Distribution illustrating the results of calculating 77-effective masses with 
formula (A. 5) for overlapped clusters from the neutral pion decays (solid lines) and the 
same data for the sample of direct photons (dotted lines), b) Distribution of 77-effective 
masses for the divided clusters from 7r decays. A mistaken contribution from the prompt 
photon due to the splitting of clusters under treatment is negligible on the level of unity 
events in the 7r -region. Smooth solid curves show the Gaussian fit with the average value 
of 133 MeV and the variation of 6.8 MeV. 
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Figure 6: Efficiency 5(7, 7) of true photon identification as a photon (solid lines), misiden- 
tification 5(7, ir°) of 7T°-meson as a photon (dotted lines), and the coefficient of the back- 
ground suppression 5(7, 7r )/e(j, 7) as a functions of E. The vector of event features is 
( J B,A 1 ,A 2 ,M 3 o,M O 4,0). 



5 Conclusion 

In this paper a neural network method is developed to separate the direct photons 
from the neutral pion background in the PHOS spectrometer of the ALICE experiment. 
The proposed algorithm is based on the analysis of the energy-profile tensor of the cluster 
calculated in its eigen vector coordinate system. The proposed method allows to construct 
effective event feature vectors consisting of a limited number of variables carrying enough 
information to train the neural network for the goals of 7/71° separation. This method 
has been applied for Monte-Carlo events in PHOS. It has been found that the probability 
of misidentification of a neutral pion as a photon is on the level of a few percent in the 
pion energy range of 3 — 120 GeV with the relatively high efficiency of the correct photon 
identification as an isolated photon in the same energy range. 
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Addendum. 

Calculation of 77-effective mass through moments M mn 

Kinematics formulae can be applied to estimate the 77-effective mass M 77 in the 
case of the divided clusters on the base of the measured energies of photons and the angle 
between their momenta. The problem is more complicated for the overlapped clusters. 
Here we follow paper [19] to express the mass M 77 in terms of the above introduced 
moments (8) under assumption that the shower profile for an isolated photon possesses 
an azimuthal symmetry with respect to the point of 7- hit. This leads to the following 
properties of the own moments M mn (i.e. moments relative to the point of 7-hit) for the 
isolated showers: M m o = Mom and M m k = for odd k. The 7-incidence inclination violates 
the azimuthal symmetry but we restore it by the mentioned compression of the cluster 
space (6). 

Due to the principle of energy additivity in the case of the overlapped showers the 
summary cluster energy density F(E 1 ,E 2 ,x 1 ,x 2 ) can be expressed in the eigen vector 
coordinate system of the cluster as follows: 



F(E 1 , E 2 , x 1 ,x 2 ) = E 1 f(x 1 - x w , x 2 ) + E 2 f(x 1 - x 20 , x 2 ) 



(A.1) 



where E 1 , E 2 are the individual energies of showers, X\ and x 2 are the coordinates of an 
arbitrary point in the cluster space, x w and x 20 are the coordinates of photon hits along 
the eigen vector e\, and f(xi,x 2 ) is the single photon shower profile. After that one can 
easily obtain from (8) 



M ln = Y J C\[E l : 

i=0 



Xiq Tfli 



+ E 2 x 



20 



rrii 



(A.2) 



where a non-zero contribution gives only terms with even i; C\ — 11/ — are bi- 
nomial coefficients, Tflifi £1X6 normalized own moments rrii n = Mj„/Moo of a single photon 
shower. Taking into account the above mentioned azimuthal symmetry of single photon 
showers, we obtain the following set of equations for the determination of energies and 
coordinates of individual showers Ei, E 2 , x w , x 20 : 



E\ + E 2 = Moo 
Ei ■ x 10 + E 2 ■ x 20 



/-I • .'7,, " E, 

E x ■ x 
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+ E 2 
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M 20 - M 02 
M 30 - 3M 12 



(A.3) 



The solution is: 



X 



M 30 - 3M 12 - M 10 (M 20 - M 02 )/M 00 



d 2 
A 
E 



(A.4) 



M 20 - M 02 - Mfo/Moo 
X 2 - 4(AM 10 - M 20 + M 02 )/M 00 
(2M 10 - XM 00 )/d 
M 00 , 

where X = x 2 o + x±o, d = x 2 o — x±o, A = E 2 — E\ and E = E 2 + E\. The variable d 
gives the estimation of distance between the hits of individual showers in the overlapped 
cluster. To calculate the effective 77-mass one should also know the distance R from the 
interaction point to the cluster center to determine the decay angle 6 between 7-quanta, 
and thus we obtain 

d 2 



M: 



77 



AE l E 2 siu 2 {e/2) = (E — A 



d? + AR 2 ' 



(A.5) 



